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Abstract 

We discuss the resummation of distributions that are singular at 
the elastic limit of partonic phase space (partonic threshold) in QCD 
hard-scattering cross sections, such as heavy quark production. We 
show how nonleading soft logarithms exponentiate in a manner that 
depends on the color structure within the underlying hard scattering. 
This result generalizes the resummation of threshold singularities for 
the Drell-Yan process, in which the hard scattering proceeds through 
color-singlet annihilation. We illustrate our results for the case of 
heavy quark production by light quark annihilation, and briefly discuss 
its extension to heavy quark production through gluon fusion. 
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1 General Formalism 

In hard scattering cross sections factorized according to perturbative QCD 
the calculable short- distance function includes distributions that are singular 
when the total invariant mass of the partons reaches the minimal value nec- 
essary to produce the observed final state. Such singular distributions can 
give substantial QCD corrections to any order in a s . 

Expressions that resum these distributions in the short-distance functions 
of Drell-Yan cross sections to arbitrary logarithmic accuracy have been known 
for some time It has also been observed that leading distributions, and 
hence leading logarithms in moment space, are the same for many hard QCD 
cross sections. This has been used as the basis for important estimates of 
heavy quark production including resummed leading, and some nonleading, 
logarithms . In this paper, we shall exhibit a method by which nonleading 
distributions may be treated, and illustrate this method in the case of heavy 
quark production through the annihilation of light quarks. 

We consider the inclusive cross section for the production of one or more 
particles, with total invariant mass Q. Examples include states produced by 
QCD, such as heavy quark pairs [[| or high-p T jets [§J, in addition to massive 
electroweak vector bosons, virtual or real, as in the Drell-Yan process. 

To be specific, we shall discuss the summation of ("plus") distributions, 
which are singular for z — 1, where 



for the production of a heavy quark pair of total invariant mass Q, with s the 
invariant mass squared of the incoming partons that initiate the hard scat- 
tering. We shall refer to z = 1 as "partonic threshold" [], or more accurately 
the "elastic limit." We assume that the cross section is defined so that there 
are no uncancelled collinear divergences in the final state. 

The main complications relative to Drell-Yan JIJ involve the exchange 
of color in the hard scattering, and the presence of final-state interactions. 
In fact, these effects only modify partonic threshold singularities at next-to- 
leading logarithm, and we give below explicit exponentiated moment-space 

1 We emphasize that by partonic threshold, we refer to cm. total energy of the incoming 
partons for a fixed final state; heavy quarks, for example, are not necessarily produced at 
rest. 
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expressions which take them into account at this level. At the next level 
of accuracy, we shall see that resummation requires ordered exponentials, in 
terms of calculable anomalous dimensions. 

The properties of QCD that make this organization possible are the fac- 
torization of soft gluons from high-energy partons in perturbation theory [||], 
and the exponentiation of soft gluon effects 0. Factorization is represented 
by Fig. 1 for the annihilation of a light quark pair to form a pair of heavy 
quarks. In this figure, momentum configurations that contribute singular 
behavior near partonic threshold are shown in a cut diagram notation |J. 
As shown, it is possible to factorize soft gluons from the "jets" of virtual and 
real particles that are on-shell and parallel to the incoming, energetic light 
quarks, as well as from the outgoing heavy quarks. Soft-gluon factorization 
from incoming light-like partons is a result of the ultrarelativistic limit [||], 
while factorization from heavy quarks, even when they are nonrelativistic, 
is familiar from heavy-quark effective theory. Once soft gluons are factored 
from them, the jets may be identified with parton distributions of the initial 
state hadrons. The hard interactions, labelled Hi and Hj in the figure, cor- 
responding to contributions from the amplitude and its complex conjugate 
respectively, are labelled by the overall color exchange in each. A general 
argument of how the exponentiation of Sudakov logarithms follows from the 
factorization of soft and hard parts and jets is given in Ref. M. 

For example, with the quark-antiquark process shown, the choice of color 
structure is simple, and may, for instance, be chosen as singlet or octet. To 
make these choices explicit, we label the colors of the incoming pair i and 
j for the quark and antiquark respectively, and of the outgoing (massive) 
pair k and I for the quark and antiquark. The hard scattering is then of the 
generic form 

H 1 = h 1 (Q 2 /fi 2 ,a s ( f i 2 ))5 ll 5 lk , (1.2) 

for singlet structure (annihilation of color) in the s-channel. For the s-channel 
octet, or more generally adjoint in color SU(N), we have, analogously 



N 2 -l 

H A = h A (Q*/^,a s (» 2 )) £ [Tf )] .. [Tf )] , (1.3) 



c=l 



with the generators in the fundamental representation. The functions 
hi are, as indicated, infrared safe, that is, free of both collinear and infrared 
divergences, even at partonic threshold. 
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Taking into account possible choices of Hj and if}, an expression that 
organizes all singular distributions for heavy quark production is 



dQ 2 dcos9* dy = $ 



:i-4) 



where y is the pair rapidity and 9* the scattering angle in the pair center 
of mass frame. The indices f and J label color tensors, such as the singlet 



1~2|) and octet Q1.3Q , into which we contract the colors of the incoming and 
outgoing partons that participate in the hard scattering. The variable S is 
the invariant mass squared of the incoming hadrons. The functions <j) a /h are 
parton densities, evaluated at scale Q 2 . The function Q contains all singular 
behavior in the threshold limit, z — > 1. 

Q depends on the scheme in which we perform factorization, the usual 
choices being MS and DIS. Note that the resummation may be carried out 
at fixed rapidity y, so long as y is not close to the edge of phase space ||. 

The color structure of the hard scattering influences contributions to non- 
leading infrared behavior. Not all soft gluons, however, are sensitive to the 
color structure of the hard scattering. Gluons that are both soft and collinear 
to the incoming partons factorize into the parton distributions of the incom- 
ing hadrons. It is at the level of nonleading, purely soft gluons with central 
rapidities that color dependence appears in the resummation of soft gluon 
effects. Each choice of color structure has, as a result, its own exponentiation 
for soft gluons M. Then, to next-to-leading-logarithm (NLL) it is possible 
to pick a color basis in which moments with respect to z exponentiate, 

n^ b J \n,y,e*,Q 2 ) = J^dzz^n^iz^e^a^Q 2 )) 

= i^%,r,QV" M * ,Q2) , (i.5) 

where the color-dependent exponents are given by 



Jo 1 — z 



dy gf ) [a s {{l-y){l-z)Q 2 )] 



o 1 - y 



+ gi ab \a s ((l-z)Q 2 )]+gi I) [a s ((l-z) 2 Q 2 ),9*] 
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+ gi J) *[a s ((l-z) 2 Q 2 ),9*} 



1.6) 



The Qi are finite functions of their arguments, and the H^ h J ^ are infrared 
safe expansions in a s (Q 2 ). g[ and are universal among hard cross 
sections and color structures for given incoming partons a and b, but depend 
on whether these partons are quarks or gluons. On the other hand, g$ 
summarizes soft logarithms that depend directly on color exchange in the 
hard scattering, and hence also on the identities and relative directions of 
the colliding partons (through 6*), both incoming and outgoing. 

Just as in the case of Drell-Yan, to reach the accuracy of NLL in the 
exponents, we need gi only to two loops, with leading logarithms coming 
entirely from its one-loop approximation, and the functions #2 and g$p only 
to a single loop. More explicitly, in the DIS scheme for incoming quarks 



MM 



with K given by 



K = C A \---^)-\n lt (1.8) 




where nj is the number of quark flavors. g<i is given for quarks by 

g^ = - 3 -C F ^. (1.9) 

L 7T 



As pointed out in [|TTJ , one-loop contributions to g% may always be absorbed 



into the one-loop contribution to g^ and the two-loop contribution to g±. Be- 
cause #3 depends upon /, however, it is advantageous to keep this nonfac- 
toring pro cess- dependence separate. We shall describe how it is determined 
below. 

First, let us sketch how these results come about [^]. After the normal 
factorization of parton distributions, soft gluons cancel in inclusive hard scat- 
tering cross sections. When restrictions are placed on soft gluon emission, 
however, finite logarithmic enhancements remain, and it is useful to separate 
soft partons from the hard scattering (which is then constrained to be fully 
virtual). Soft gluons may be factored from the hard scattering into a set of 
Wilson lines, or ordered exponentials, from which collinear singularities in 
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the initial state are eliminated, either by explicit subtractions or by a suitable 
choice of gauge ||. Assuming that the lowest-order process is two-to-two, 
there will be two incoming and two outgoing Wilson lines.Q The result, illus- 
trated in Fig. lb, is of the form, H^Su, summed over the same color basis 
as in eq. (|1.4| ) above. 

The resulting hard scattering and soft-gluon functions both require renor- 
malization, which is organized by going to a basis in the space of color ex- 
changes between the Wilson lines. The renormalization is carried out by a 
counterterm matrix in this space of color tensors. For incoming and outgoing 
lines of equal masses, such analyses have been carried out to one loop in M, 
and ]13| and to two loops in a related process in fTI |. For an underlying 



partonic process a + b — > c + d, we then construct an anomalous dimension 
matrix p^ - *^ where the indices I and J vary over the various color ex- 
changes possible in the partonic process. The soft function SVj then satisfies 
the renormalization group equation [RJ 



(4 +/3(9) l) s " 



Svj> . (1.10) 



This resummation of soft logarithms is analogous to singlet evolution in 
deeply inelastic scattering, which involves the mixing of operators, and hence 
of parton distributions. The general solution, even in moment space, is given 
in terms of ordered exponentials which, however, may be diagonalized at 
leading logarithm. For the resummation of soft logarithms in QCD cross 
sections, the same general pattern holds, with mixing between hard color 
tensors. Leading soft logarithms, however, are next-to-leading overall in mo- 
ment space, which allows the exponentiation (|1.5|) at this level. 

Of course, the analysis is simplest for external quarks, and most compli- 
cated for external gluons. It is also possible to imagine a similar analysis 
when there are more than two partons in the final state. This would be 
necessary if we were to treat threshold corrections to pp — > QQ + jet, for 
instance, but we have not attempted to explore such processes in detail. 

Given a choice of incoming and outgoing partons, next-to-leading loga- 
rithms in the moment variable n exponentiate as in (|1.5| ) in the color tensor 
basis that diagonalizes Y^fj^ cd \ with eigenf unctions A/. The resulting soft 

2 In Drell-Yan and other electroweak annihilation processes, there is a pair of incoming 
lines only. 
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function g$ is then simply 

#[a s ((l - z) 2 Q 2 ), 8*} = -\j[a s ((l - z) 2 Q 2 ), 9*} , (1.11) 

where the eigenfunctions are complex in general, and depend on the directions 
of the incoming and outgoing partons, as shown. 



2 Applications to qq — >■ QQ 

These considerations may be illustrated by heavy quark production through 
light quark annihilation, 



Q(Pa) + q{Vb) -> Qipi) + Q(P2) ■ 



Following Ref. [[OJ, we define invariants 

*i = ip a -Pif ~ m 2 , u 1 = (p b -p 1 ) 2 - 
with m the heavy quark mass, which satisfy 

s + ti + u\ = 



m 



S = {Pa+Pbf 



(2.1) 



(2.2) 



(2.3) 



at partonic threshold. In this case, as in elastic scattering || [14|], the anoma- 
lous dimension matrix is only two-dimensional. 

In a color tensor basis of singlet and octet exchange in the s channel, the 
anomalous dimension of eq. ( l.lUj ) is, 



r 2 i 

Tl 2 
^22 



ry 

S -C F (Lp+l+7Ti) 
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(2.4) 
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where subscript 1 denotes singlet, and 2 octet. Lp is the familiar velocity- 
dependent eikonal function 

1 - 2m 2 Is (l-3\ j , 

with p = J 1 — 4m 2 /s. The matrix depends, as expected, on the directions 
of the Wilson lines, which may be reexpressed in terms of ratios of kinematic 
invariants for the partonic scattering. For arbitrary f3 and fixed scattering 
angle, we must solve for the relevant diagonal basis of color structure, and de- 
termine the eigenvalues. However, V is already diagonalized in the s-channel 
octet-singlet basis at "absolute" threshold, (3 — 0, and, more generally, when 
the parton-parton cm. scattering angle is 9* = 90° (where U\ — tx). 

It is of interest, of course, to compare the one- loop expansion of our 
results to known one-loop calculations, at the level of next-to-leading order. 
We may give our result as a function of z, since the inverse transforms are 
trivial. They are found in terms of the Born cross section, the one-loop 
factoring contributions of g± and and 1^2 • In the DIS scheme the 

result is 



^ J) (^Wi,ti,s) (1) = CTBorn ~~ ~ ~ I C[ 

ij n I z I 



21n(l-z) + ! 



8]u ( j- J -2-2L^ + 21nM 



2 



+c A 



(2.6) 



Here [i is the factorization scale, and the logarithm of s/^j 2 describes the 
evolution of the parton distributions. This result cannot be compared directly 
to the exact one-loop results of 0] for arbitrary j3, where the singular behavior 
is given in terms of the variable S4, with 

s 4 = (p 2 + k) 2 -m 2 ^2p 2 -k, (2.7) 

where k = p a + — pi — p 2 is the momentum carried away by the gluon. At 
partonic threshold, both s 4 and (1 — z) vanish, but even for small s 4 , angular 
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integrals over the gluon momentum with S4 held fixed are rather different 
than those with 1 — z ~ 2(pi + p^) ■ k/s held fixed. 

Nevertheless, the cross sections become identical in the (3 —> limit, 
where we may make a direct comparison. Near s = 4m 2 , we may identify 
2m 2 (l - z) = s 4 , and (2.6) reduces to the (3 -> limit of eq. (30) of @. 
It is also worth noting that even for (3 > 0, the two cross sections remain 
remarkably close, differening only at first nonleading logarithm in the abelian 
(Cfi) term, due to the interplay of angular integrals with leading singularities 
for the differing treatments of phase space. 

As for the Drell-Yan cross section, our analysis applies not only to abso- 
lute threshold for the production of the heavy quarks (/3 = 0), but also to 
partonic threshold for the production of moving heavy quarks. When j3 nears 
unity, however, the anomalous dimensions themselves develop (collinear) sin- 
gularities associated with the fragmentation of the heavy quarks, which in 
principle may be factored into nonperturbative fragmentation functions. 

In summary, we have illustrated the application of a general method for 
resumming next-to- leading logarithms at partonic threshold in QCD cross 
sections. Possible extensions include, of course, heavy-quark production 
through gluon fusion, and dijet production. Extensions to multijet produc- 
tion are also possible. We reserve estimates of the phenomenological impor- 
tance of these nonleading terms to future work, but we hope that whether 
they give small contributions or large, the method will improve the reliability 
of perturbative QCD calculations for hard scattering cross sections. 
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Figure Caption 

Cut diagram illustrating momentum configurations that give rise to thresh- 
old enhancements in heavy quark production, (a) General factorization the- 
orem. Away from partonic threshold all singularities in the "short-distance" 
subdiagram H/S cancel; (b) Expanded view of H/S near threshold, showing 
the factorization of soft gluons onto eikonal (Wilson) lines from incoming and 
outgoing partons in the hard subprocess. Hj and Hj represent the remaining, 
truly short-distance, hard scattering. 
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